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Abstract 



The digital filter has an interpolated signal value (SW) provided as the sum of a given number of weighted 
sample values, with each weighting factor determined by an interpolation function with a corresponding 
number of linear sections. The interpolation function may have a succession of linear sections of constant 
gradient, each extending over a time interval of half the time spacing between the successive sample values 
which are weighted by the weighting factors. 
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Digitale Filterstruktur 

v Urn die Abtastrate um ein beliebiges Verhaltnis zu yer- 
andern, wird bei eiriem digitalen Filter zur zeitlichen Inter- 
polation diskreter Signale eine Funktion zur Gewichtung 
der vorgegebenen Stutzwerte benutzt, die abschnittswei- 
se linear definiert ist. Das digitale Filter zeichnet sich 
durch etnen geringen Aufwand, insbesondere bei der Pro- 
gram mie rung der Funktion, auf. Die Anmeldungsgebiete 
liegen bei der Verarbeitung von digitalen Audio- und Vi- 
deosignalen. 
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Beschreibung * 

Die Erfindung bezieht sich auf ein digitales Filter zur zeitlichen Interpolation diskreter Signale. Solche Signale erhalt 
man beispielsweise durch Abtastung analoger Audio- oder Videosignale. 
5 Bei der digitalen Signalverarbeitung werden neben den durch Abtastung gewonnenen Signale haufig zusatzliche Si- 
gnale benotigt, die das urspriingliche analoge Signal charakterisieren. Zusatzliche Signale lassen sich aus den Abtastwer- 
ten durch Interpolationsverfahren rnit Interpolationsfiltern berechnen. Im allgemeinen Fall ermittelt ein Interpolations fil- 
ter aus n Eingangswerten m Ausgangswerte. m kann dabei groBer oder auch kleiner als n sein. 1st die Anzahl der Aus- 
gangswerte kleiner als die Anzahl der Eingangswene, liegt eine Dezimation vor. Bei einer Dezimation urn den Faktor 2, 
' 10 also wehn gilt m = 0,5 • n wird von den.n Eingangswerten nur jeder zweite Wert ausgegeben. Die Errnitdung von Aus- 
gangswerten ist komplizierter, wenn beispielsweise aus funf Eingangswerten vier Ausgangswerte erzeugt werden sollen. 
In der Regel sind Ausgangswerte erwiinscht, die voneinander jeweils den selben konstanten zeitlichen Abstand aufwei- 
sen. Solche aquidistanten Werte erhalt man nicht, wenn einfach jeder funfte Wert nicht an den Ausgang weitergeleitet 
wird. Vielrnehr sind die Zeitpunkte der vier Ausgangswerte gegenuber den Zeitpunkten der fiinf Eingangs werte verscho- 
15 ben. Im Gegensatz dazu fallen die Ausgangswerte bei einer Dezimation \im einen ganzzahligen Faktor mit den entspre- 
chenden Eingangswerten zusammen. 
I Ein Verfahren zur beliebigen Anderung von Abtastraten ist in dem Artikel "Digital Methods for Conversion Between 
Arbitrary Sampling Frequencies" von Tor A. Ramstad in IEEE Transactions on Acoustics, Speech, and Signal Proces r 
sing, vol. ASSP-32, Nr. 3, Juni 1984, Seiten 577 bis 591 beschrieben. Die zugehorigen Schaitungen werden dort als hy- 
20 bride Systeme bezeichnet; Sie bestehen aus einem digitalen Vorfilter und einem nachgeschalteten analogen Interpolator 
(Res ampler). Aufgrund der groBen Zahl arithmetischer Operationen, insbesondere Multiplikatiohen, ist eine schaltungs- 
technische Realisierung sehr aufwendig! '- 

Aufgabe der Erfindung ist es daher, ein digitales Filter zur zeitlichen Interpolation diskreter Signale anzugeben, das bei 
vergleichbaren Interpolationsergebnissen mit einem geringeren Schaltungs auf wand als herkoinmliche Filter realisierbar 
. ™ 25 ist; Gelost wird diese Aufgabe mit eihem digitalen Filter mit den Merkxnalen des Anspruchs 1 , .< , . 

Die Erfindung hat* deh ydrieii, dafi die mterpolationsf unktion abschhitts weise linear definiert ist;, Jeder Abschnitt wird 
also durch : einen einf achen mamemadscheriZusammenhang beschrieben. 5 ■ ' . > . 

Die Implementierung der Interpolationsfunktion rnit einer Programmiersprache ist besonderes einfach; da die Funk- 
tion in Abschnitte unterteilt ist. Die Symmetrie-Eigenschaften der Interpolationsfunktion kbnnen in vorteilhafter Weise 
30" bei der Prograrrimierung beriicksichtigt werden. Einander entsprechender Abschnitte unterscheiden sich nur im \brzei- 
chen ihrer Steigung. .* • . t , 

Mit dem Parameter k lassen sich Breite und Dampfung der Storbereiche im Amplitudehgang des Filters, die' aufgrund 
der Periodizitat des Filters beim Ubergang einer Periode zur folgenden auftreten, einstellen. Bei Werten fur k < 1 ist die 
Dampfung schwacher und das Filterverhalten breitbandiger als bei k = 1. 
35 Die Eigenschaften des digitalen Filters konnen,weiter verbessert werden, wenn seine Irnpulsantwort mit einer Recht- 
. eckfunktionpassende hange gefaltet wird. 

Vorzugsweise wird das digitele Filter mit einem linearphasigen Interpolatibnsfilter kombiniert. Das Interpolationsfilter 
beispielsweise verdoppelt die Abtastrate und stellt dem digitalen Filter so wohl die urspriinglich vorgegebenen Werte als ■ 
, auch die durch Interpolation gewonnenen Zwischen werte zur Verfuguhg. Im digitalen Filter wird dann die Abtastrate urn 
40 ein beliebiges Verhaltnis verandert. Durch geeighete Wahl der Konstanten der Ubertragungsfunktion des Interpolations- v 
filters kann eine optimale Anpas sung an das digitale Filter erreicht werden. 

Weitere vorteilhafte Aus- und Weiterbildungen sind in Unteranspruchen gekennzeichnet. . * 
Die Erfindung wird hachfolgend anhand von Figuren riaher erlautert/ Es zeigen; 
Fig.* 1 die Ubertragungsfunktion der Interpolationsfunktion, f -, / ■ . 

45 Fig. 2 das zeitliche Verhalten-der Interpoladpnsmnktion in bezug;auf die.Stutzwerte, ' • ' \\ . , ( 

'"*"* - Kg.. 3 eine FUterstruktur mit dem digitalen Filter und 

r ^ jJB3g._4 e in linear phasig es Interpolationsfilter ' - ^- I ^ _ ■ •_ ■ 

Das erfindungsgeniaBe digitale Filter berechnet aus vorgegebenen Stutzwerten, die beispielsweise durch die Abta- 
stung eines analogen Signals bekannt sind ? einen Signal wert SW zu einem beliebigen Zeitpunkt zwischen zwei benach- 
5b barren Stutzwerten so, daB er dem wirkliehen Wert des Signals zu diesem Zeitpunkt m6glichst gut entspricht. Fiir diese 
Interpolation werden Stutzwerte verwendet, die sich in zeitlicher Nachbarschaft zu dem zu interpolierenden Signalwert 
SWbefinden/ 

Jeder der verwendeten Stutzwerte wird mit einem Faktor bewertet, der sich aus einer Interpolationsfunktion a(t) ergibt. 
Bei jedem zu interpolierenden Signalwert S W wird die Interpolationsfunktion a(t) so gewahlt, daB sie beziiglich der Zeit 
, 55 zu einem Zeitpunkt ti, zu dem der Signalwert SW zu berechnen ist, symmetrisch ist. Jeweils zu den Zeitpunkten der 
Stutzwerte laBt sich ein entsprechender Cjewichtungsfaktor ermitteln, mit dem der jeweilige Stutz wert multipliziert wird. 
Die Summe der Produkte aus diesen Multiplikationen ergibt den Signalwert SW zu dem Zeitpunkt ti. 

GemaB Fig. 1 ist die Interpblationsfunktion a(t) abschnittsweise linear definiert. Sie weist zwischen 0 und 2,5 Zeitein- 
heiten funf Abschnitte Bl, B2, B3, B4 und B5 auf. Jeder Abschnitt ist 0,5 Zeiteinheiten lang. Die ersten beiden Abschnitt 
60 Bl, B2 sind monoton steigend, die beiden letzten Abschnitte B4, B5 sind monoton fallend. Uber dem Abschnitt B3 ist 
die Funktion konstant. Die Amplitude A ubersteigt zu keiner Zeit t den Wert 1 . 

Filter weisen ublicherweise im Amplitudengang eine Periodizitat von 2T auf. Beim Ubergang einer Periode zur fol- 
genden treten Storbereiche auf, also bei 2T, 4T, 6T Mit einem Parameter k der Interpblationsfunktion a(t), der Werte 

zwischen einschlieBlich -1/2 und +1/2 annehrnen kann, konnen Dampfung und Breite der Storbereiche eingestellt wer- 
65 den. Bei einem Wert von k < 1 ist die Dampfung schwacher und das Filterverhalten breitbandiger, also von einem Wert 
von k = 1 . In Fig. 1 ist der Parameter k der Interpolationsfunktion a(t) jeweils fur einen Wert von k eingetragen, der =1, < 
1 oder > 1 ist. Da die Interpolationsfunktion a(t) abschnittsweise definiert ist, lassen sich die einzelnen Abschnitte jeweils 
getrennt beschreiben. Im einzelnen gilt dabei: 
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l/4-k/8 + (2-k)dT/4,wennfurtgilt: to + 0<t < Iq+1/2 
1/2 + k • dT/2, wenn fiirt gilt: to + 1/2 < t < k) + 1 
1/2 + k/4, wenn fur t gilt: (<, + 1 < t < to + 1,5 
1/2 - k ■ dT/2, wenn fur t gilt: tb + 1,5 < t < to + 2 

1/4 - k/8 - (2 - k) dT/4, wenn fur t gilt: to + 2 < t < t 0 + 2,5. 5 

AUe Zeitangaben beziehen sich auf eine beliebig definierte Zeiteinheit. Bei der obigen Beschreibung gibt to ein Start- 
zeitpunkt vor, der 1,25 Zeiteinheiten vor dem Zeitpunkt tj des zu interpolierenden Signalwerts SW liegt. dT ist die zeit- 
liche Verschiebung des zu interpolierenden Signal wertes SW zu seinem nachsten benachbarten vorgegebenen Stutzwert 

Die Linie der Kurve der Interpolationsfunktion a(t) ist fur den Parameter k = 1 durchgezogen, fur den Parameter k < 1 
mit langen Strichen gestrichelt und fur k > 1 mit kleinen Strichen gestrichelt. 10 

Bei der Darstellung nach Fig. 1 wurde to = 0 gewahlt. D.h. ein zu interpolierenden Signalwert SW wiirde beim Zeit- 
punkt ti = 1,25 Zeiteinheiten liegen. 

Mit Hilfe einer Sprungfunktion s(t - 1) lafit sich die Interpolationsfunktion a(t) auch in geschlossener Sehreibweise an- 
geben. Die Sprungfunktion s(t - 1) beschreibt einen Einheitssprung, der bis zum Zeitpunkt t = 1 den Wert 0 und ab die- 
sem Zeitpunkt fur alle folgenden Zeitpunkte den Wert 1 aufweist. Die Interpolationsfunktion a(t) lautet dann: 15 
" a(t) = (t - kt/2) • s(t) + (-2 + 3k) • (-1 + 2t)/4 • s(t - 1/2) + (k - k • t) • s(t - 1) + k • (3/2 - 1) • s(t - 3/2) + ((-2 + 3 • k) • (-2 
+ t))/2 • s(t - 2) - ((-2 -f k) • (5 + 2 • t) /4 ■ s(t - 2,5). 

Wie bereits erlautert bestimmt der.Parameter k, der Werte zwischen einschlieBlich -1/2 und +1/2 annehmen kann, die 
Breite und Dampfung der ungeraden Storbereiche. 

Wie kurz vprher erlautert, werden die vorgegebenen Stutzwerte mit der Interpolationsfunktion a(t) gewichtet Gra- 20 
phisch betrachtet muB die Interpolationsfunktion a(t) auf der Zeitachse so verschoben werden, daB eine gedachte'Sym- 
metrieachse der Interpolationsfunktion durch den Zeitpunkt tj, zu dem der zu interpolierende Signalwert SW zu bestim- 
men ist, 1 auf t. Auf grand des angenahert trapezformigen Verlauf der Kurve der Interpolationsfunktion a(t) stent jede ge- 
dachte Symmetrieachse senkrecht zur Zeitachse t. In Fig. 2 sind die Zeitpunkte t k _i, t k und t^+i yon Stutzwerten n k _x, rit 
und n k+ i.. eingetragen. Der zu dem Zeitpunkt ii des zu interpolierenden Signalwerts SW zeitlich nachst gelegene Stutzwert 25 
ist n k zu dem Zeitpunkt i k . E)er zeitliche Abstahd zwischen diesen beiden Zeitpunkten wird als Verschiebung; dT bezeich- . 
net. Diese Verschiebung dT 1st bei der abschnittsweisen Beschreibung der Interpolationsfunktion a(t) die ^Variable. Zu 
den Zeitpunkten der Stutzwerte ri k _i, n k und n k+ i werden die entsprechenden Werte der Interpolationsfunktion a(t) ermit- 
telt. Beim Stutzwert n k _i erhalt man den Wert a(t k _ x ), beim Stutzwert n k den Wert aft^) und beim Stutzwert n k+1 den Wert 
a(t k+ i ) . Der zeitliche Ab stand zwischen zwei unmittelb ar aufeinander folgenden S tutzwerten betragt eine Zeiteinheit. Die 30 
Betrage dieser Stutzwerte werden mit dem entsprechenden Wert a(t) multipliziert. Die Summe meser drei Produkte gibt 
ein Betrag, der dem gesuchten Signalwert SW zu dem Zeitpunkt t^ entspricht. Es ist zu beachten, daB n k _ x , n k und n k+1 je- 
weils eine Benennung des entsprechenden Stiitzwertes und nicht die jeweiligen Amplituden werte bezeichnen. 

Das digitale Filter berechnet also zu einem beliebigen Zeitpunkt ti'aus Stutzwerten, die zu diesem Zeitpunkt Benach- 
bart liegen, den Signalwert SW und gibt diesen aiis. I . 35 

Fig. 3 zeigt eine bevorzugte Kombination von Filteranordnungen unter Verwendung des erfindungsgemaBen digitalen 
Filters, daB hier mit F2 bezeichnet ist. Signale werden mit einer Frequenz %r abgetastet und einem linearphasigen Inter- 
polationsfilter Fl zugefuhrt. Dort wifd die Abtastrate verdoppelt, d. h. zu jedem urspriinglichen Abtastwert wird'durch 
Interpolation ein Zwischenwert erzeugt. Die urspriinglichen Abtastwerte und die im Interpolationsfilter berechneten 
Zwischenwerte, werden an das digitale Filter F2,weitergegeben.-Dort kann die Abtastrate urn ein beliebiges Verhaltnis 40 
verandert werden. Die so ermittelteri Signalwerte werden von einem Dezimationsfilter F3 gefiltert, umeine Bandbegren- 
zung zu erzielen. Dies ist dann notwendig, wenn die Abtastrate der Signalwerte groBer ist als die Abtastrate f^. An ei- 
nem Ausgang des Dezimationsfilters F3 stehen die ermittelten Werte mit der Abtastfrequenz foirr zur weiteren Signal- 
verarbeitung zur Verfugung. 

Fig. 4 zeigt ein Ausfuhrungsbeispiel fur das Jinearphasige Interpolationsfilter. Es handelt sich dabei unrein Filter 5. 45 
Drdnung mit 5 Verzogerungsgliedem z" 1 . Ein Eingangssignal x durchlauft eine bestimmte Anzahl von Verzogerungsglie- 
dern bevbr es mit sich selbst urn eine bestirnmte,Zeit verzogert addiert wird. Die Summe der Verzogerungen der beiden 
Summanden entspricht stets 5 Verzogerungsgliedem z" 1 . Auf diese Weise erhalt man 3 Zwischenwerte, die mit Faktoren 
a, b und c multipliziert werden. Aus der Summe dieser Produkte erhalt man einen Ausgangswert Yl , der dem zuvor ge- 
' nannten Zwischenwert entspricht. Durch eine bloBe Verzogerung des Eingangswertes X bekommt man einen weiteren 50 
, Ausgangswert Y2. Die Verzogerung des weiteren Ausgangswertes Y2 gegenuber dem Eingangswert X ist erforderlich, 
da der Ausgangswert Yl durch die Interpolation einer zeitlichen Verzogerung unterliegt und die zeitliche Reiherifolge 
des Zwischenwerte zu den urspriinglichen Werten eingehalten werden muB. 

Das linearphasige Interpolationsfilter. Fl nach Fig. 4 laBt sich beschreiben durch: 

• 55 

H(Z) = a - b • Z" 1 + c • Z 2 + c • Z" 3 - b • Z^ + a • Z" 5 . 

Um eine gute Anpassung an das nachgeschaltete digitale Filter F2 zu erreichen, gelten fur die Faktoren a, b und c fol- 
gende Bedingungen: 

0 < a < 1/4 60 
0 < b < 1/2 
c = 1/2 + b - a. 

Eine besonders gute Anpassung erhalt man bei der Wahl von: 
a = 1/32, b = 5/32 und c = 5/a. - 

Gute Ergebnisse erzielt man auch mit: 65 
a = 1/32, b = 9/64 und c = 39/64. 

Selbstverstandiich konnen auch linearphasige Interpolationsfilter hohere Ordnung mit dem digitalen Filter F2 kombi- 
niert werden. 
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Patentanspriiche 

1. Digitales Filter zur zeitlichen Interpolation diskreter Signale, dadurch gekennzeichnet, daB das Filter einen zu 
interpolierenden Signalwert (SW) als Surnme aus mit jeweils einem Faktor (a(t)) gewichteten vorgegebenen Stiitz- 
werten (n k _!,.n k , n k+1 ) berechnet, die zu dem Signalwert (SW) benachbaft sind, wobei jeder Faktor (a(t)) durch eine 
Interpolationsfunktion der Zeit, die uber die Zeit abschnittsweise linear definiert ist, zu einem Zeitpunkt (tt-i, t k , 
tk+i) des entsprechenden Stiitzwertes (n^, n k , n n+1 ) bestimmt ist. , 

2. Digitales Filter nach Anspruch 1 , dadurch gekennzeichnet, daB die Interpolationsfunktion jeweils in einem Zeit- 
abschnitt (bl, b2 v . . b5) linear mit konstanter Steigung ist, der halb so lang ist wie der zeitliche Abstand zwischen 
zwei unrnittelbar aufeinander folgenden Stutzwerten (n*, n k+ i). 

3. Digitales Filter nach Anspruch 1 oder 2, dadurch gekennzeichnet, daB die Interpolationsfunktion auf der Zeit 
achse (t) beziiglich des Zeitpunktes (tp des zu interpolierenden Signalwertes (SW) symmetrisch ist. 

4. Digitales Filter nach einem der Anspriiche 1 bis 3, dadurch gekennzeichnet, daB die Interpolationsfunktion ein 
Zeitintervall (B3) mit einem konstanten Funktionswert aufweist. 

5. Digitales Filter nach einem der Anspriiche 1 bis 4, dadurch gekennzeichnet, daB die Interpolationsfunktion be- 
schreibbar ist- durch: 

a(t) = (t - kt/2) ■ s(t) + (-2 + 3k) (-1 + 2t)/4 ■ s(t - 1/2) + (k - k • t) ■ s(t - 1) + k ■ (3/2 - t) • s(t - 3/2) + ((-2 + 
3-k).(-2 + t))/2-s(t-2)-((-2 + k),(5 + 2t));/4-s(t^ . 

wobei a(t) der Funktionswert zu einem Zeitpunkt t ist und 1,25 Zeiteinheiten vor den Zeitpunkt des Signalwertes 
(SW)giit:t = 0, 

k ein von der Zeit. unabhangiger Parameter ist und . . 

s(t - 1) eiri Einheitssprung ist, dej bis zum Zeitpunkt t =1 den Wen 0 und ab diesjsm Zeitpunkt den Wert 1 aufweist 
% , Digitales Filter nach einem der Anspriiche 1 bis 4, dadurch gekennzeichnet^ daB die InterpolaUonsfunktion zu 
Zeitpunkten t beschreibbar ist durch , '' » .-•> '" 

■ 1/4 ^ + (2 ^k) 4 dTTi, wehn ftir t gUt: 1^ + 0 < it <' to 4- 1/2 . v : ; ^ ' : ' * ."' ■/'• : . V\. 
1/2 + k v dT/2,' wenn fur/£jlt: .to;.+ 1/2'- < v t*"< -•t 0 , .+'l , " 4 ■/ 

1/2 + k/4, wenn fiir t gilt: to + 1 < t < tb+1,5 . *- • . .* . '" 

1/2 - k ■ dT/2, wenn fiir t gilt: to + 1,5 < t < to + V ■ - . : -\ 

1/4 - k/8 - (2 - k) dT/4, wenn fur t gilt: ^ + 2 < t ^ to + 2,5, 

wobei der zeitliche Abstand zwischen zwei unrnittelbar aufeinariderfolgenden vorgegebenen Stutzwerten (n k ; n k+ i) 
eine Zeiteinheit betragt, ' : , . ' "• 

' to ein Stanzeitpunkt ist, der 1,25 Zeiteinheiten vor dem Zeitpunkt (tj) des zu interpolierenden Signalwertes (SW) 
liegtund . / , - * ' 

dT die zeitliche Verschiebung des zii interpolierenden Signalwertes (SW) zu seihem nachsten benachbarten vorge-. 
gebenen Stiitzwert (nk) darstellt. 

7 . Digitales Filter nach einem der Anspruche 1 bis 6, dadurch gekennzeichnet, daB ein linearphasiges Interpolati- 
onsfilter (Fl) einen Teil der vorgegeben Stiitzwerte (n^) liefert. 

8. Digitales Filter nach Anspruch 7, dadurch gekennzeichnet, daB die Ubertragungsfunktion des Interpolationsfil- 
ters (Fl) beschreibbar ist.durch: . , , 

H(Z) = a-tVZ- 1 + c , - 

wobei a, b und c Konstanten sind. , L ■ 1 - 

9., Digitales Filter nach Anspruch 8, dadurch gekennzeichnet, daB fiir die Konstanten gilt: > 

o<a<i/4 •-• * ; . . '\ " r " • ■■ \ /■ 

c= 1/2 + b-a. 
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